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Abstract - We study the unitary time evolution of the order parameter of a quantum system after 
a sudden quench in the parameter which drives the transition. By mapping the dynamics onto 
the imaginary time path integral in a film geometry we derive the full mean-field non-equilibrium 
phase diagram for a one-component order parameter. The recently discovered non-equilibrium 
transition is identified with the shifted critical point in films and therefore it is generally expected 
to occur in more than one spatial dimension. We also find that anharmonic oscillations of the 
order parameter are a general feature of the mean-field dynamics after the quench. 



Introduction. — The experimental realization [T] of 
many-body cold atomic systems evolving in the absence 
of dissipation triggered a significant theoretical activity 
aimed at understanding the fate of a quantum interacting 
system initially prepared in a state which is not a Hamilto- 
nian eigenstate. A global quantum quench is the prototyp- 
ical instance: the system is prepared at time i = in the 
ground state of a Hamiltonian Hq and then it is let evolve 
with a Hamiltonian H which differs from Hq everywhere 
in space (e.g., due to a sudden quench of an external field). 
The absence of a solid and general theoretical framework 
for studying the ensuing non-equilibrium evolution moti- 
vated the investigation of simple one-dimensional models 
with exactly solvable dynamics (see, e.g., Refs. [SHS], but 
this list is far from being exhaustive). Several aspects of 
their evolution have been explored both analytically and 
numerically [7], as reviewed in Ref. [8]. Only few studies 
of models in spatial dimension d > 1 are presently avail- 
able ^EM- 

One of the central issues concerning the time evolution 
after a quench is whether a stationary state exists at long 
times. It is commonly believed [11I31[SJ[71[TT1[T2] that, in 
case it does, it is characterized by correlation functions 
which resemble those of either a thermal state or a so- 
called generalized Gibbs ensemble. However, a station- 
ary asymptotic state does not always exist since the order 
parameter of the system might display persistent anhar- 
monic oscillations, as it occurs in the case of the mean-field 
(MF) approximation for superfluid bosons [T3HT5] , for the 



Hubbard model |16| and for some exactly solvable pairing 
models ^7M\- 

A general approach to the study of quench problems, 
proposed in Ref. [3] , consists in mapping the real-time evo- 
lution of a d-dimensional quantum system onto a bound- 
ary problem in d -I- 1 dimensions in imaginary time. With 
this framework the quantum evolution governed by critical 
Hamiltonians [3|[20] or by massive integrable field theories 
in d = 1 [2Tj, has been successfully investigated. 

In this paper we consider such a mapping for a quench 
within a generic effective Landau-Ginzburg (LG) Hamilto- 
nian, which is expected to capture the universal properties 
of a rather large class of quantum critical points, encom- 
passing some of the cases mentioned above. For a one- 
component order parameter, i.e., within the Ising univer- 
sality class, we present a remarkably simple complete MF 
solution for the order parameter evolution which, as we 
shall see, accounts comprehensively for all the qualitative 
features observed in real-time MF analysct0- The main re- 
sult of our investigation is that the novel " dynamical tran- 
sition" (which has no equilibrium counterpart) reported in 
Refs. [T4l[T6] is the manifestation, in real time, of a rather 
well-known effect in confined critical systems, i.e., the crit- 
ical point shift in films |22j . Remarkably, the universality 
of this phenomenon allows us to conclude that the dynam- 



^ In the case of systems with n-componcnt vector order parameter 
and 0{n) symmetry — such as those discussed in Refs. I13ltl7l — 
the analysis reported here applies to the only component of the order 
parameter which docs not vanish when the O (n) symmetry is broken 
by the boundaries. 
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ical transition is not a MF artifact but it occurs even in 
the presence of fluctuations, an important issue which was 
left open by all the analytical approaches considered so 
far in the literature. In addition, wc find that anharmonic 
oscillations are general features of the MF evolution and 
we characterize quantitatively their features (form, am- 
plitude, and period) in terms of the quench parameters. 
However, the existence of the transition beyond MF does 
not guarantee that these oscillations do actually withstand 
the effects of quantum fluctuations. Experiments and ad- 
ditional theoretical studies should be able to provide a 
conclusive answer to this subtle and difficult question. 

The mapping to the film geometry. Following 
Ref. [3], consider a d-dimensional system prepared in a 
state \^o) which is not an eigenstate of the Hamiltonian 
H controlling the time evolution. The expectation value 
of a local operator 0{x) at time t is given by the path- 
integral 



iHt-eH\ 



(1) 



H 



where Z = (v^ole^ l^'o) and the damping factors e 
have been introduced in order to make the path-integral 
representation of {0{t,x)) absolutely convergent. Equa- 
tion ([ij can be represented by an analytically continued 
path integral in imaginary time r in which the field ^' 
takes the boundary values for t = ri.2 = ±e — it. 
The operator O is inserted at t = 0. The width of the 
film is T2 ~ Ti = 2e. In the course of the calculation ti^2 
must be considered as real numbers and only at the very 
end they are continued to their effective complex values. 
In this way, the real-time non-equilibrium evolution of a 
d-dimensional system is mapped onto the thermodynam- 
ics of a d -I- 1-dimensional field theory in a film geometry 
in which the initial state jvPo) plays the role of a bound- 
ary condition at both boundaries. This approach can be 
used to study the time-evolution of a generic operator O. 
However, in what follows wc will concentrate on the order 
parameter (j)(t,x). 

When the system is close to a (quantum) phase transi- 
tion, we can take advantage of the renormalization-group 
(RG) theory of boundary critical phenomena [23| in order 
to describe the system in terms of the LG Hamiltonian for 
the order parameter (/)(t, x) in a film |23] : 



d'^.T / dr 

Tl 



+ / A'^x [B{4>{t^,x))+B{^(t2,x))] , (2) 



[d = (V,9i-)] where 6(0) = \c(j)^ — hs(j> represents the 
boundary density, r controls the distance from the critical 
point, and g > ensures the stability for r < 0. c is 
related to the surface enhancement, i.e., to the difference 
between the interaction strengths within the boundaries 
and within the bulk, whereas is a symmetry-breaking 
field acting at both confining surfaces. The values of c and 





Fig. 1: Order parameter profiles close to the boundary, for 
hs ^ and c > (left), c < (right). The RG invariant 
boundary conditions, i.e., (t>s — and = ±cxd, are effectively 
imposed at t — ti — |A|. For hs large enough the profile looks 
like the right one. 



hs determine the actual boundary conditions for the field 
(j). This is easily understood in the semi-infinite geometry 
T2 — oo, taking into account the MF boundary condition 
-~dr<t>s +c(f)s — hs~0 implied by Eq. ^ for (ps ^ '/'(t'i, x), 
whereas (j)f, = (j){T — > oo,a;) attains its bulk value i/if, = 
\/~-^rJg for r < 0, and 0b = for r > 0. For hg = and 
(ps 7^ 0, 101 increases upon leaving the boundary if c > 0, 
as shown by the left panel of Fig. [TJ The order parameter 
profile extrapolated outside the film vanishes for r = ri — A 
(see Fig. [1]), where A is the so-called extrapolation length 
A ~ 1/c for c large enough [23]. If c < 0, instead, |0| 
decreases upon leaving the surface, as shown by the right 
panel of Fig. [U and the (negative) extrapolation length A 
determines the position r = ri + A at which the MF order 
parameter profile diverges. In order to have 0s 7^ with 
hs — 0, one needs either r > and — c > r^/^ or r < 0, 
which correspond to 0h = or 0b 7^ 0, respectively. A 
large enough hs causes an enhancement of |0| close to 
the surface, so that 0^ does not vanish independently of 
the value of (pi, and the profile looks like the one depicted 
on the right panel of Fig. [TJ 

Accounting for fiuctuations requires a RG analysis [23j 
which yields the following conclusions for translational in- 
variant boundary condition: for hs = 0, the possible fixed- 
point values of c are c* = 00, 0, — cxj, with c* = ±00 sta- 
ble and c* = unstable. In particular c* = +00 leads 
to Dirichlet (D) boundary condition for at the bound- 
aries of the film, i.e., 0s flows imdcr RG towards 0* = 0. 
The fixed-point with c* ~ —00 describes both the cases 
(— c > r*, /is = 0) [53] and (c < 00, hs ^ 0) and leads to a 
divergence of the order parameter profile upon approach- 
ing the smface: 0(t — >■ 0,2;) ~ t~^/'^ (with (3 and v the 
standard critical exponents), i.e., 0s flows to 0* — ±cxj, 
referred to as ± boundary conditions. The kind of asymp- 
totic boundary conditions can be also expressed in terms 
of 0s and 0h. If 0s = 0, the corresponding boundary con- 
dition is Dirichlet. On the other hand, if 0s 7^ 0, the 
order parameter profile looks like on the left or right of 
Fig. [TJ depending on whether 0s < 4>b or (ps > 4>b, which 
correspond respectively to D and + boundary condition. 

Going back to the two boundaries problem, the most im- 
portant feature of the phase diagram is the critical point 
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shift [22] ■ This is easily understood in the case of clas- 
sical critical points for which r = T — Tc, being Tc the 
critical temperature in the bulk. For DD boundary condi- 
tions, the suppression of order at the surfaces prevents the 
film from ordering and therefore (f) = everywhere even 
slightly below T^. However, upon further decreasing T, the 
increasingly strong tendency towards ordering overcomes 
the effects of the boundaries and the whole film orders. 
This occurs at the shifted critical point T = Tc{L) < Tc, 
which depends on the thickness L of the film, and is such 
that Tc{L) - Tc ^ L-^l'' . For boundary conditions, 
such a point is additionally shifted off-coexistence and it 
can be reached by applying a non-zero bulk magnetic field 
which overcomes the order at the surface. 

Mean-field order-parameter profiles - Let us consider 
the LG Hamiltonian (21) in a film geometry of thickness 
L = T2— Ti and introduce the rescaled MF order-parameter 

profile 0(r) = T\f^'^^^^ ~ which satisfies the 

equation —"0" -I- fV" + V''^ = 0, where f = rL^ . For DD 
boundary conditions, ■i/'(0) = V'(l) — 0, whereas for ++ 
boundary conditions, — > 0) ~ and il)(x — > 1^) ~ 
(1 — x)~^. The differential equation can be solved ana- 
lytically. For H — h boundary conditions and f > fc = — tt^ 
(high-temperature phase H) one finds [25] 



il){x) 



: , . ^ fT^r. dn(2Kx]k) 



(3) 



where K = K{k) is the complete elliptic integral of the 
first kind and the elliptic modulus k = k{f) solves f = 
(2A')^(2fc^ — 1). dn and sn arc the standard Jacobi elliptic 
functions. For f < fc (low-temperature phase L), instead, 
one finds [25l 



ip++{x) = 2V2K [sn{2Kx; fc)]" 



(4) 



where fc satisfies f = — (2/v)^(fc^ + 1). 

For DD boundary conditions, ij.' vanishes for f > fc- For 
f < fc, instead, the film orders and the order parameter 
is dS] 

iPooix) ^2V2K ksn{2Kx;k) , (5) 

where now fc solves f = —{2K)'^{k'^ + 1). In Fig. Hj^a) we 
plot for the cases of +-|- and DD boundary conditions. 

The analytic continuation. For the purpose of 
extracting the asymptotic behavior of (f)[t) for t ^ tq we 
replace |^o) by the appropriate RG-invariant boundary 
condition Iv^q) to which it flows. The difference is taken 
into account, to leading order, by assuming that the RG- 
invariant boundary conditions are not imposed at r = ti,2 
but at r = Ti^2 T tq, where tq is the absolute value of 
the extrapolation length. In the quantum non-equilibrium 
problem, tq is of the order of the correlation length in the 
initial state [3] . The effect of introducing tq is simply to 
replace e by e -I- tq and the limit e — > 0+ can now be taken, 
resulting in a film of effective width L = 2tq. In order 
to determine the time evolution of the order parameter 
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Fig. 2: (a) Order parameter profiles ^p{x) across a film with DD 
and H — h boundary conditions. The profiles denoted by L and 
H correspond, respectively, to the low- and high-temperature 
phases, (b) Value of the order parameter in the middle of 
the film as functions of f. L and H indicate the values of f 
corresponding to ip{x) in panel (a). In both panels, the dotted 
lines refer to the bulk order parameter. 



(j> one has simply to continue analytically the expressions 
provided above for -0 to the complex position t = tq -\- it 
within the film. Such a continuation is straightforwardly 
done by using standard formulas for the Jacobi functions 
[24]. For a quench with H — h boundary conditions to a 
point below rc = —tt'^/{4:Tq), one finds 



[L] 



r/\rc\ 



a Kdn{Kt/TQ;k'), 
-(2i^/7r)2(fc2 + l) < 



-I, 



(6) 



where fc' = Vl — fc^ and a = 2\/3/ (To^/g)- 4>+{t) oscillates 
between aK and akK with period Tfc = 2tqK' /K, where 
K' = K{k'). For ++ boundary conditions, but with r > 
rc- 



[H] 



-1 



(7) 



(j)+{t) ^ ak' Kcn{Kt/To;k') , 
r/\rc\ = (2if/^)2(2fc2 - I) > 

(j)+{t) oscillates with period 2Tfc between —ak'K and 
ak' K. The qualitative change in the evolution below and 
above the boundary point occurs at r = rc via an expo- 
nential relaxation 



(f>+{t) cx [cosh(7ri/2ro)]" 



(8) 



which is similar to the result at criticality in one dimen- 
sio and is characterized by a decay time Toxp = 2tq/'k. 
Upon approaching rc. 



Tk{r^rc 



,ln(24/|I-r/rc|) 



(9) 



and therefore the periods Tk and 2Tk of the oscillations 
([6]) and ([7]) — respectively above and below rc — diverge 



^We mention that from the mathematical standpoint, anhar- 
monic oscillations in one spatial dimension are absent because for 
critical phenomena in a two-dimensional film geometry the Jacobi 
functions which characterize the mean-field expressions encountered 
so far for the order parameter profile are replaced by trigonometric 
functions (see, e.g., 1271 ). The analytic continuation of the latter 
involve hyperbolic functions, which result in an exponential relax- 
ation towards the corresponding asymptotic values. Physically, this 
is a consequence of the absence of ordered phases in classical one- 
dimensional systems and therefore the two low-temperature regions 
disappear. 
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Fig. 3: Left: Inverse period of the oscillations of in unit 

of T,7^ as a function of r/|rc|. For r — >■ Tc the period diverges 
logarithmically. Right: Maximum and minimum value of 0+ (t) 
(up to the multiplicative constant a) as a function of r/\rc\. 

logarithmically with a prefactor set by the exponential 
decay time Tgxp. For DD boundary conditions we have 
instead 

[H] (j}D{t)^Q for r>re, (10) 

and 




so that (/>d(0 oscillates with period Tk between aK and 
akK. For r < r^, 4>D{t) differs from (j)+{t) only by a half- 
period shift. Figure [3] shows the period and the amplitude 
of the oscillations of 0+ whereas in Fig. U] we present the 
three typical non- vanishing evolutions. In the case of ++ 
boundary conditions a mean-field analysis of the critical 
evolution (r = 0) was presented in Ref. [3], where anhar- 
monic oscillations were reported. Here we have extended 
the analysis to the entire phase diagram with generic r 
and boundary conditions, with a careful characterization 
of the period of the emerging anharmonic oscillations as 
a function of the quench parameters. This careful study 
allowed us to highlight the appearance of the dynamical 
mean-field transition and reveal its very nature. 

Discussion. — Our results depend explicitly on the 
parameter tq, which determines the value of Tc cx Tq'^ . In 
Ref. [3] it has been argued that tq is proportional to the 
spatial correlation length in the initial state, the inverse 
of the initial mass gap. Within the mean-field LG formal- 
ism ^0 is proportional to l^ol"^/^, where rg is the initial 
value of r in Eq. ([2]). Accordingly, ttq is proportional to 
r/lrol and therefore Tc (x tq, but with unknown propor- 
tionality factors. Although we cannot fix the exact value 
of Tq, the non- equilibrium phase diagram, which gives the 
appropriate asymptotic boundary conditions correspond- 
ing to a certain choice of the quench parameters rg and r, 
can be understood on physical grounds. In the case of dis- 
ordered initial state rg > 0, one always expects Dirichlct 
boundary conditions in the film, independently of the final 
value r. If such a value exceeds the boundary transition 
point Tc, the order parameter </> vanishes at all times. Gon- 
versely, for r < rc, one has anharmonic oscillations the fea- 
tures of which depend on rg through rg oc |ro|~^/^. When 




Fig. 4: Order parameter (in arbitrary units) as a function 
of t/ro: (a) (j)+{t) for r < is always positive, (b) (f)+if) for 
r > rc oscillates symmetrically around 0. (c) 4'D{i) for r < r^. 
Note that (poit) = for r > r^. 



the quantum system is in an initial state with rg < 0, the 
asymptotic boundary conditions depend on the value of 
r. From the general features of the profiles in the film 
(sec Fig. (2), we expect that if 0g > (p^q, the corresponding 
asymptotic boundary conditions is ++, whereas it is DD 
for (j>Q < (/)oq, where (j>cq and (po correspond, respectively, 
to (j)b and (j}s- The crossover between these two regimes 
occurs on the no-quench line r = rg. All these considera- 
tions lead to the non-equilibrium phase diagram sketched 
in Fig. [5] (where, in the spirit of our analysis, r, rg, and Tg 
are considered as independent parameters). Remarkably, 
this phase diagram accounts for all the qualitative features 
found via real-time MF approaches [T4 l ll6 j . In particular, 
in the Bosc-Hubbard model [Tl|, for quenches within the 
ordered (superfluid) phase, ((/"^(i)) displays two types of 
oscillations separated by the point of "dynamical tran- 
sition" at which an exponential relaxation occurs, with 
decay time Tgxp- These two different kind of oscillations 
actually correspond to " phases" (II) and (III) in Fig. [51 as 
it can be realized by comparing the time evolution of the 
square modulus of the order parameter presented in Fig. 2 
(A) and (G) of Ref. [14] (dashed lines) with the one which 
can be readily inferred for cfP from Fig. |4] (a) and (b), re- 
spectively. Upon approaching the dynamical transition, 
the period Tj- of these oscillation — expressed in terms of 
the corresponding Texp — was found [14] to diverge accord- 
ing to the leading term of Eq. ^ . In addition, it turned 
out that for a fixed initial value Uq of the coupling constant 
U of the Bose- Hubbard model the dynamical transition oc- 
curs at a value U'^{Uo) = (C/q -I- Uc)/2 of the coupling U 
after the quench [14], where Uc is the value of U at which 
the model undergoes a Mott-superfluid transition in equi- 
librium. (The notation here is slightly different from the 
one used in Ref. [E]). In terms of the distance r oc U —Uc 
from the critical point, the dynamical transition therefore 
occurs at rc oc U'^{Uo) — Uc ^ {Uq — Uc)/2 oc rg where 
rg oc [/g — Uc is thc value of r right before the quench. 
The fact that rc oc rg provides an additional evidence of 
the correctness of the picture presented here, which al- 
lows us to conclude that the point of dynamical transition 
highlighted in Ref. [13] is nothing but r = r^ i.e., the 
shifted critical temperature in the corresponding classical 
film. This implies that the dynamical transition is not a 
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Fig. 5: Mean-field non-equilibrium phase diagram. We identify 
four "phases": (1) r > Tc and ro > 0, where (jf>(t) vanishes iden- 
tically; (II) r > Tc and ro < 0, where (j}{t) shows anharmonic 
oscillations around zero; (III) and (IV) where (j}{t) shows an- 
harmonic oscillations with a non-zero average. These latter two 
regions are separated by the no-quench line r = ro and, qualita- 
tively, the corresponding (j){t) differ only by a half-period shift. 
Upon approaching the line r — Vc, the period of oscillation di- 
verges according to Eq. ^ and (j}{t) decays exponentially to 
zero with decay time Tcxp. 



MF artifact and therefore it has to be expected generically 
for d > 1. This could also be the case for the dynam- 
ical transition points with anharmonic oscillations found 
in the MF evolution of the Hubbard [T5l[T6] and pairing 
models [l7l[T8]. first, the period of the anharmonic oscil- 
lations of the order parameter diverges upon approaching 
it and, second, it separates a region [(I) + (II)] within 
which the asymptotic time average of the order parame- 
ter vanishes from a region [(III) -I- (IV)] within which it 
does not. Fluctuations beyond MF are likely to affect and 
completely wipe out these oscillations, resulting either in 
a relaxation of the order parameter towards some asymp- 
totic value or even to some chaotic dynamics in which 
oscillations (damped or not) with various frequencies are 
involved, as suggested by the results of Ref. [14] . Within 
the present approach, however, it is difficult to establish 
whether these fluctuations are relevant only below an up- 
per critical dimensionality (as it is the case of equilibrium 
critical phenomena) or, instead, they destroy the oscilla- 
tions in generic dimensionality. Nonetheless, the RG and 
scaling analysis of the classical system in the film geome- 
try — which accounts for critical fluctuations — suggests 
that the presence of the dynamical transition is robust in 
any dimensionality (7^ 1) and that this transition, as its 
classical counterpart, is characterized by the second of the 
two features mentioned above. However, the critical-point 
shift Tc is affected by fluctuations and indeed it scales as 
~ Tq where v is the equilibrium critical exponent 
governing the divergence of the order parameter correla- 
tion length ^ of the d+ 1-dimensional classical system with- 



out boundaries. This fact provides an additional example 
of how the dynamical behavior after a quantum 
quench may display some of the universal features of the 
corresponding classical equilibrium model. 
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